Quantum Computing by Two-Dimensional NMR using Spin- and Transition-Selective 

Pulses 
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Quantum computing using two-dimensional NMR has recently been described using scalar cou- 
pling evolution technique [J. Chem. Phys., 109, 10603 (1998)]. In the present paper, we describe 
two-dimensional NMR quantum computing with the help of selective pulses. A number of logic gates 
are implemented using two and three qubits with one extra observer spin. Some many-in-one gates 
(or Portmanteau gates) are implemented. Toffoli gate (or AND/NAND gate) and OR/NOR gates 
are implemented on three qubits. Deutsch-Jozsa quantum algorithm for one and two qubits, using 
one extra work qubit, has also been implemented using selective pulses after creating a coherent 
superposition state, in the two-dimensional methodology. 

Key Words: Two-dimensional NMR; quantum computation; logic gate; Toffoli gate; Deutsch 
Jozsa algorithm. 



INTRODUCTION 



Quantum computing exploits the intrinsic quantum 
nature of physical systems and is therefore more powerful 
than classical computing for a certain class of problems. 
While some problems like factorization gain exponential 
speed up (1), some others like database-search gain poly- 
nomial speed up (2). A quantum computer works on 
two-level quantum systems known as 'quantum bits' or 
'qubits'. The fact that qubits can exist in a coherent 
superposition of basis states is exploited in a quantum 
computer. However, retaining such a coherent superpo- 
sition for long enough time is one of the major hurdles 
in quantum computation. Since nuclear spins in liquids 
are efficiently isolated from the rest of their environment, 
the coherent superposition lasts for several hundreds of 
milliseconds. Therefore, high resolution NMR of weakly 
coupled nuclear spins in liquids is one of the potential 
candidates for realizing such a quantum computer. Nu- 
clear spins with 1=1/2 have generally been used since 
spins with J>l/2 relax much faster. Quantum comput- 
ing using one-dimensional NMR has been demonstrated 
by various workers (3-12). 

Recently, quantum computing has been demonstrated 
using two-dimensional NMR (13). In this method, the 
states of the 'computation' spins (or input spins) are 
encoded by the transitions of an extra 'observer' spin 
(Fig.l). The observer spin is first allowed to evolve for 
a time tx during which, the input qubits remain in their 
initial state (Fig. 2). After the frequency labelling period 
tx, the computation is performed on the input qubits. 
The observer spin is again allowed to evolve for a time t2 



and detected. The experiment starts from a mixed input 
state, the computation is performed simultaneously on all 
the input states, and one gets a mixed output state. One 
can also create a superposition of input qubits in the be- 
ginning of the experiment by applying a it/ 2 pulse on the 
input qubits and subsequently killing the coherence by a 
gradient pulse (13). Various steps in two-dimensional 
NMR namely, preparation, evolution, mixing and detec- 
tion have a close correspondence with the steps in quan- 
tum computing namely, creation of initial states, labeling 
of the initial states, computation and reading of output 
states, respectively (13) (Fig. 2). In this respect, two- 
dimensional NMR appears to be the 'method of choice' 
for quantum computing. 

An advantage of the two-dimensional method is that 
it resolves all the input and the output states and corre- 
lates them. This correlation between input and output 
states in the two-dimensional experiment makes the re- 
sult of the computation graphic. For example, a SWAP 
gate which exchanges the states of two qubits, can be im- 
plemented in NMR by selectively interchanging popula- 
tions of zero-quantum levels. Since the one-dimensional 
NMR spectrum of homonuclear spins after the execu- 
tion of SWAP gate is indistinguishable from that of an 
equilibrium spectrum, the one-dimensional SWAP gate 
was earlier demonstrated after the creation of a non- 
equilibrium state (12, 14). However, as will be shown 
here (Fig. 3), and has been demonstrated earlier (13), 
the two-dimensional method starts from the equilibrium 
state and yields a spectrum characteristic of the SWAP 
gate. 

Two-dimensional quantum computation has been car- 
ried out so far by using scalar coupling evolution tech- 
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nique (13). This paper describes a two dimensional NMR 
implementation of several two and three qubit gates 
with one extra observer spin, using spin- and transition- 
selective pulses. The first implementation of the Deutsch- 
Jozsa quantum algorithm using two-dimensional NMR is 
also demonstrated here. One of the advantages of the 
selective pulse method is that it makes the computation 
simple and straightforward. 

EXPERIMENTAL 
A. Logic Gates 

A weakly coupled 3-spin system (I0I1I2) can be viewed 
as a combination of two 2-spin subsystems with the ob- 
server spin (I ) being in state or 1 (Fig.l). 



001 _.„--- 




Figure 1. Energy level diagram of a 3-spin system (I0I1I2) 
showing two 2-spin subsystems corresponding to states and 
1 of observer spin Io . The transitions of Io are labelled by the 
states of input qubits (I1I2). 

There are four observer spin transitions labelled as 1 1 , 
10, 01 and 00, which correspond to the states of the in- 
put qubits (I1I2) in these transitions. During the compu- 
tation (mixing) period, various transitions of the input 
qubits are inverted yielding various gates as listed in Ta- 
ble 1. The results of 24 one-to-one reversible gates for the 
3-spin system having one observer and two input qubits 
are shown in Fig. 3. The NOP gate is implemented by do- 
ing No OPeration during the computation period. In this 
gate, each input state corresponds to same output state 
after the computation. Various NOT gates are imple- 
mented by inverting one or both of input qubits (Ii T2) , 
using spin-selective n pulses. XOR (or control-NOT) and 



XNOR gates are implemented by inverting two similar 
transitions of the same input spin, one in each subsys- 
tem. We label XOR and XNOR gates as XOR1, XNOR1 
and XOR2, XNOR2 depending on whether the result of 
the operation is stored on the spin Ii or I 2 respectively. 
For example, XOR1 gate is implemented by inverting 
transitions 001-011 and 101-111 of Ii and the result of 
the gate is stored on Ii (Table 1, Fig.3). A SWAP gate 
can be implemented by selectively interchanging the pop- 
ulations of zero quantum levels of each subsystem (Table 
1, Fig.3). This can be achieved by the use of a cascade 
of three transition-selective, non-commuting tt pulses on 
regressively connected transitions during the computa- 
tion period (12, 14). It may be noted that while all the 
gates of ref.(13) are reproduced here, with several new 
gates added, the methodology used for the computation 
is different except for the NOP and NOT gates. 

Similar to classical Boolean expressions, quantum 
gates can be reduced to a sequence of simple unitary op- 
erations. Although a few basic gates such as NOT, AND, 
OR, XOR arc sufficient to carryout a given computation, 
other gates are also useful since they may reduce the 
number of pulses. Therefore 'portmanteau gates' which 
perform more than one operation have been implemented 
in one-dimensional NMR using selective pulses (12). 
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Figure 2. Pulse scheme for the two-dimensional NMR im- 
plementation of quantum computing. The close correspon- 
dence between two-dimensional NMR and quantum comput- 
ing is shown. I is the observer spin and Ii, I2 are the input 
qubits. The gradient pulse G z eliminates the unwanted trans- 
verse magnetization before the computation. 
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Figure 3. X H NMR spectrum (b) of 2,3-dibromo-propionic acid (a) in CDCI3 recorded on a Bruker DRX-500 spectrometer 
at 300K. (c) shows observer spin (Io) spectra corresponding to various gates implemented using spin- and transition-selective 
pulses. Pulse scheme used is shown in Fig. 2 and various transitions of input qubits (Ii, I2) inverted during the computation are 
listed in Table 1. The spin-selective pulses were 10ms long and the transition-selective pulses were 100-300ms long. The phase 
of the computation pulses were cycled through (x, —x) to suppress the distortions due to pulse imperfections. All experiments 
were carried out in the time domain with 256 ti values and 256 complex data points along t2 and with 2 scans for each ti 
point. Zero filling to 512x512 complex data points was done prior to the 2D Fourier transformation. All plots are shown in 
magnitude mode. 
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Various logic gates with Boolean algebra, truth table and operations performed 
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a \s,t) represents the state of input qubits (Ii,I 2 ). 

b Order of transitions are important as all transitions may not commute. 



Similar gates have been implemented here using two- 
dimensional NMR (Table 1, Fig. 3). Interchanging pop- 
ulations of double quantum levels using a cascade of 
non-commuting ir pulses on progressively connected tran- 
sitions leads to SWAP+NOT gate (Table 1, Fig.3). 
SWAP+XOR and SWAP+XNOR gates have been im- 
plemented by inverting two pairs of regressively con- 
nected transitions in each case (Table 1, Fig.3). Sim- 



ilarly, SWAP+NOT+XOR and SWAP+NOT +XNOR 
gates have been implemented by inverting two pairs of 
progressively connected transitions (Table 1, Fig.3). The 
last six gates in Fig.3 are direct combinations of two 
gates. It may be noted that the 24 gates shown in Fig.3 
form a complete set of 2-qubit one-to-one mappings. 
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Figure 4. 19 F and : H NMR spectra (b) of 2,3-difluro-6- 
nitrophenol (a) in CDCI3 (with one drop of D2O to induce 
the exchange of the hydroxy proton and hence to suppress 
its coupling to fiourine nuclei) recorded on a Bruker DRX- 
500 spectrometer at 300K. (c) shows 19 F spectra of observer 
spin Io corresponding to various gates. Pulse scheme used 
is shown in Fig. 2 and various transitions of control qubit 
Ii inverted during the computation period are described in 
the text. The spin-selective pulses were 1ms long and the 
transition-selective pulses were 200ms long. The phase of the 
computation pulses were cycled through (x, —x) to suppress 
distortions due to pulse imperfections. All experiments were 
carried out in the time domain with 1024 ti values and 256 
complex data points along I2 and with 2 scans for each ti 
point. Zero filling to 1024x1024 complex data points was 
done prior to the 2D Fourier transformation. All plots are 
shown in magnitude mode. 



Fig. 4 shows several 3-qubit gates implemented on a 4- 
spin system, using selective pulses. Once again, No Oper- 
ation during computation period yields NOP gate and in- 
verting spin Ii yields the NOT(Ii) gate. The more inter- 
esting ones are - Toffoli gate (or AND/NAND gate) and 
OR/NOR gate. The operations of Toffoli and OR/NOR 
gates (13,15) are respectively, 



\s, t,u) — > \s © (t A u),t, u) 
and 

\s,t,u) — > |s © (t V u), t, u), 



(1) 



(2) 



where © = addition modulo 2, A = AND, V = OR and 
s, t, u are the states of the control spin Ii and the input 
spins I2 and I3. The Toffoli gate is a universal gate for 
reversible computation. This gate can be implemented 
by inverting the control spin transitions 011 and 111, by 
using a pair of transition-selective 7r pulses. Inverting all 
transitions of the control spin except the transitions 010 
and 110 leads to OR/NOR gate (Fig.4). 



B. Deutsch-Jozsa algorithm 

The Deutsch-Jozsa (DJ) algorithm demonstrates the 
power of quantum computing over classical computing 
(16). In this algorithm, we consider functions from N-bit 
domain space to 1-bit range space. A function is called 
constant if it gives the same output for any input, and is 
called balanced if it gives one output for half the number 
of inputs and another for the remaining half. Given an N- 
bit binary function which is either constant or balanced, 
classically 2^ N ^ 1 > + 1 operations are required to determine 
whether the function is constant or balanced, whereas DJ 
algorithm requires only a single evaluation. The DJ algo- 
rithm has been implemented using one-dimensional NMR 
by several research groups (10-13). 

In the Cleve version of the DJ algorithm (17), a bi- 
nary function / is encoded in a unitary transformation 
by the propagator U/ by including an extra work qubit 
such that 



\r)\s)—*\r < 



'f(s))\s), 



(3) 



where |r) and \s) are respectively the states of work qubit 
(Io) and input qubit (Ii). The four possible functions for 
the single-bit DJ algorithm are listed in Table 2. 

TABLE 2 
The four possible binary functions 
(/1-/4) for 1-bit DJ algorithm 
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Figure 5. Pulse scheme for the two-dimensional NMR im- 
plementation of DJ algorithm. Io is the work qubit and Ii,±2 
are the input qubits. 



This algorithm has been implemented on the thermal 
equilibrium state and does not require the creation of 
a pure initial state (10,12). The two-dimensional pulse 
scheme used for implementing DJ algorithm is shown in 
Fig. 5. The experiment begins with both qubits in a su- 
perposition of states, achieved by a non-selective (Tr/2) y 
pulse. This is followed by an evolution period ti, prop- 
agator U/ and detection period t2. The transformations 
corresponding to /i and fa are respectively, a unity op- 
eration and a spin selective n x pulse on the work qubit. 
The transformations corresponding to fa and fa are im- 
plemented by transition-selective ir x pulses respectively 
on the 10-11 and 00-01 transitions of the work qubit (12). 
The results of the algorithm for all the 4 functions are 
shown in Fig. 6. The constant or balanced nature of the 
function is identified by the presence or absence of signal 
from the input qubit (Ii) (Fig. 6). The expected signals 
are also shown schematically in Fig. 6. For two of the 
functions fa and fa, the expected signals are calculated 
in Table 3. The spin-states (|r)|s)) before computation 
can be paired (connected by a curved line) in such a way 
that each pair represents a transition of the input qubit 
(Ii). Each transition of the input qubit is labelled by 
the state of the remaining spin i.e., work qubit. Cor- 
responding output pairs of spin-states after computation 
can be calculated by using expression [3] . From each out- 
put pair, the corresponding transition of the input spin 
after computation can be identified. If two spins have 
flipped in the output pair, then the transition becomes 
non-observable and will be labelled as zero quantum (ZQ) 
or double quantum (DQ). 
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Figure 6. 1 H NMR spectrum (a) of 5-nitro furaldehyde (b) 
in CsD6 on a Bruker DRX-500 spectrometer at 300K. The 
results of DJ algorithm for various functions /1-/4 in Table 
2 are shown in (c). I is the work qubit and Ii is the input 
qubit. Expansions of the Ii parts of the spectra are shown 
in dotted boxes. The expected pattern is also shown in each 
case. The expected patterns for fa and fa are described in 
Table 3. The spin-selective pulses were 10ms long and the 
transition-selective pulses were 100ms long. The phase of the 
computation pulses were cycled through (x, —x) to suppress 
the distortions due to pulse imperfections. All experiments 
were carried out in the time domain with 256 ti values and 
512 complex data points along t2 and with 2 scans for each ti 
point. Zero filling to 512x512 complex data points was done 
prior to the 2D Fourier transformation. All plots are shown 
in magnitude mode. 
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TABLE 3 

The input-output correlations for the functions 
/ 2 and / 4 of 1-bit DJ algorithm. 
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Implementing DJ-algorithm on two input qubits re- 
quires three qubits including one work qubit. The algo- 
rithm can be described as 



\r)\a)\t)%®f(a,t))\s)\t), 



(4) 



where \r), \s) and \t) are the states of the work qubit (Io) 
and of two input (Ii,l2) qubits respectively. The eight 
possible two-bit binary functions are listed in Table 4. 

TABLE 4. 

The eight possible binary functions (/i-/s) 
for the 2-bit DJ algorithm 
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The pulse scheme is same as in Fig. 5. Once again, the 
transformations corresponding to /i and fi are respec- 
tively unity operation and a spin selective ir x pulse on the 
work qubit. The unitary transformations encoding the 
six balanced functions /V/s are implemented by selec- 
tive pulses on the transitions of the work qubit, taken two 
at a time i.e., [0, 0, n, tt], [tt, tt, 0, 0], [tt, 0, ir, 0], [0, tt, 0, tt], 
[7r,0,0,7r], and [0, 7r,7r, 0] where denotes no pulse on 
that particular transition (10,12). The results of the al- 
gorithm for all the 8 functions are shown in Fig. 7. 
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Figure 7. The results of DJ algorithm on 2,3-dibromo- 
propionic acid (Fig.3a,b) in CDCI3 for various functions fi- 
fs listed in Table 4. Io is the work qubit and I1J2 are the 
input qubits. Only expansions of the Ii and I2 parts of the 
spectra are shown. Expected pattern is also shown in each 
case for comparision. The expected patterns for fe and fi 
are described in Table 5. All experiments were carried out on 
a Bruker DRX-500 spectrometer at 300K. The spin-selective 
pulses were 10ms long and the transition-selective pulses were 
100-300ms long. The phase of the computation pulses were 
cycled through (x, —x) to suppress the distortions due to 
pulse imperfections. 1024 complex data points along t2 were 
acquired for 256 ti values with 2 scans for each ti point. Zero 
filling to 1024x1024 complex data points was done prior to 
the 2D Fourier transformation. All plots are shown in mag- 
nitude mode. 



Once again, a function is constant only if signals from 
all the input qubits are present, otherwise the function is 
balanced (Fig. 7). Table 5 describes the input-output cor- 
relation for /2 and fi in a 2-qubit D J algorithm. Here the 
s- and t- transitions are labelled by the states \r)\t) and 
\r)\s) respectively. Other details of Table 5 are similar to 
that of Table 3. 
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TABLE 5. and Technology, New Delhi are gratefully acknowledged. 

The input-output correlations for the functions We also thank Prof. K. V. Ramanathan and Dr. G. A. 
f 2 and f 4 of 2-bit DJ algorithm Naganagowda of SIF for discussions. 
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CONCLUSIONS 

Using two-dimensional NMR quantum computing 
method, several 2- and 3-qubit gates have been imple- 
mented on 3 and 4 weakly coupled spin systems and 1- 
and 2-qubit DJ algorithm on 2 and 3 weakly coupled spin 
systems by utilizing spin- and transition-selective pulses. 
The use of selective pulses over scalar evolution has sev- 
eral advantages as well as some disadvantages. The ad- 
vantages of selective pulses are that they lead to simple 
logic, need only a few pulses and work very well. The dis- 
advantages are that, one needs long low power r.f. pulses 
during which relaxation and r.f. inhomogeneity effects 
degrade the performance of the selective pulses and ideal 
selectivity may be difficult to achieve. Furthermore, for 
the selective pulse experiments, one needs resolved tran- 
sitions. However, many of these difficulties can be over- 
come by the use of heteronuclear spins such as 15 N, 13 C, 
19 F, 31 P and 1 H, where the magnitudes of the spin-spin 
couplings are large, the transitions are well spread out 
and selectivity is easier to achieve using pulses of shorter 
durations. Further improvement can be achieved by ori- 
enting molecules in liquid crystal matrices. Attempts are 
continuing in this direction in our as well as in other lab- 
oratories (18,19). 
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